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1. Introduction 

Inflationary paradigm has dominated cosmology since long ago and recently its 
predictions have attained excellent observational evidence. However at present there 
exists no preferred concrete inflationary scenario based on a convincing realistic high 
energy physics model. It is therefore interesting to explore every single theoretical 
possibility which leads to inflation along with specific predictions. 

In the standard chaotic inflation scenario [I] the massive inflaton field rolls down 
its potential. If its initial amplitude is much larger than Planck mass Mpi, the Hubble 
friction is large and inflaton field is in the "slow-roll" regime where its potential energy 
dominates over its kinetic energy. In this case the energy density during "slow-roll" is 
nearly a constant and the equation of state is nearly that of the cosmological constant. 
The exit from such inflationary regime occurs when the Hubble parameter becomes 
of order of inflaton mass. When this happens the inflaton field starts to behave as 
ordinary matter and later decays into light particles. This simple model describes 
current observations rather well. There are numerous examples which develop this 
idea further j2j. The predictions of all these models are quite similar, so that even the 
recent accurate measurements [3J did not put too much pressure on the viability of most 
of these models. Probably the most serious constraint is the presently unobserved and, 
thus, necessarily small fraction of isocurvature fluctuations in the primordial spectrum. 
These fluctuations generically arise in multi-field models. The other constraint comes 
from the nongaussianity in the primordial spectrum |3j. Future experiments on the 
detection of the gravitational background radiation [Sj left after inflation may impose 
further significant constraints on the inflationary models [6j. 

Most of the inflationary models can be classified as potential induced, i.e. inflation 
is driven by the potential terms in the Lagrangian of a scalar field. This is, however, not 
the only theoretical approach to implement the inflation. There is a distinct possibility 
that the inflationary regime was induced by the kinetic energy. The kinetic terms, of 
course, have to be of a nonstandard form, since for models with standard kinetic terms 
and without potential the universe is always decelerating. It is, however, possible to 
implement an inflationary evolution of the early Universe if one assumes a nonstandard 
Lagrangian. A novel way to generate an inflationary regime using the above approach, 
has been proposed in and dubbed " k -inflation" . The Lagrangian in this model was 
taken in the form of expansion in orders of the standard kinetic term X = \ (<9 M 0) 2 : 

£ = K(4>)X + L(4>)X 2 + ... =p{<t>,X). (1) 

The system described by the Lagrangian of this form was shown to dynamically enter the 
inflationary regime starting from rather generic initial conditions. There is an analog 
of a "slow-roll" regime, where the systems has an adiabatically changing, nearly de 
Sitter equation of state. The model possesses the energy-momentum tensor of the 
same structure as for the perfect fluid with the energy density e and the pressure p. 
Therefore the hydro-dynamical intuition is useful for the analysis. At the end of k- 
inflation the equation of state w = p/e changes dynamically from w ~ — 1 to the 
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"ultra-hard" equation of state w = 1 . Thus, exit from k-inflation occurs in a natural 
graceful manner. 

In this paper we propose a new model for inflation where the Lagrangian depends 
only on the second derivatives of a scalar field. The simplest scalar that can be built 
from a scalar field and its second covariant derivatives is the D'Alembert operator 
acting on 0. We will restrict ourselves to this case. In fact, a large class of Lagrangians 
of the form 

C=Q(Of>) 1 (2) 

where Q(D0) is a convex function, will be shown to drive an inflationary evolution of 
the universe from rather generic initial conditions. This theory has to be viewed as 
an effective field theory where the first derivative terms are suppressed. One can, for 
example, suppose that in the Lagrangian |IJ) the functions K(<p) and L((j)) vanish in some 
limiting cases. Then the higher order derivative terms may dominate the dynamics of 
the system for some period of its evolution. The Lagrangians of this type were also 
considered in jHHEl in various contexts^. 

We will dub this model analogously to k-inflation as "Box-inflation" or for 
simplicity " B -inflation" . The system described by this "truncated" Lagrangian does 
not have an exit from inflation, so after discussing the properties of this system during 
the inflationary stage we propose simple modifications of the model which provide a 
graceful exit from B-inflation. One natural way to implement this exit is to revive the 
kinetic terms which depend on and X. For the proposed modifications, the system 
will generically exit to the state with w ~ 0. This is similar to the graceful exit in 
the case of chaotic inflation. The perturbation theory for the B-inflation as well as a 
detailed analysis of the observational consequences will be developed elsewhere [T2] . 



2. The Model 



We will first consider the action containing the D'Alembert operator □ = g^V^V 1 " 
(where V is the covariant derivative) acting on a scalar field minimally coupled to 
gravity: 

2~ 



S 



Sg + 



d X\f—g 



-M 



+ M 2 M 2 Pl Q 



(3) 



M 2 M Pl 

Here Q is an arbitrary convex function (Q" > const > 0, later we will see that this 
inequality is important). We will only consider the case when Q(0) = in order not to 
have the cosmological constant being inserted by hand. This underlines the fact that the 
inflationary regime in our model is a result of dynamics of the system. In addition it is 
natural to assume that the system is symmetric with respect to the Z 2 transformations 
— > —(f). Thus we can take only even functions Q into consideration. The scale M, at 

% Note that, despite of similarities with the model from Ref. ^Hji our system does not have any 
potential and cannot be reduced to the noninteracting two field model pi] as it was done in the 
reference above. Thus the physical content of our model is completely different from that of Ref. [10| . 
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this point, is an arbitrary parameter. For the purposes of this paper its value is not very 
important. We will leave the question of the value of this scale for a future discussion. 
The "strange" choice of combination of masses for the Q-term is due to the following 
reason. Writing the action S<f, as a Taylor series in (□</>) we find: 

"Q"(o)^ )2 + Q (4) (o) ^ 4 



Sa 



d Xy/—g 



□c 



MJ,,M 6 



+ 



(4) 



2M 2 v Ty ' 4! 

Now we can see that the initial choice of combination of masses leads to a simple form 
of the first non-zero term in Taylor expansion of action S^. Let us mention that, in 
general, we will not rely on the expansion (J3J). There is another reason for choosing the 
combination M 2 M V \ in Q. Since the action is invariant under the "shift" symmetry 
— ► + 0o, where 0o is an arbitrary constant, the energy density is insensitive to the 
value of 0. It is, however, sensitive to the values of derivatives of 0. Thus it is natural 
to make the D-operator dimensionless using one scale, whereas the value of the field 
is made dimensionless with some other scale. The scale M can be viewed as the 
physical cutoff scale in our model. The other scale with which is made dimensionless 
we are free to choose at our convenience. We will choose this scale to be a Planck scale. 
With this choice, it follows from the Friedmann equation that the Hubble parameter is 
naturally measured in units of M. 

Let us return to Eq. Q and perform the following substitution: 

-> M P ,0, x^ -> M~V. (5) 

Note that the space-time coordinates and the Hubble parameter H are made 
dimensionless with the same scale. The operator □ is also dimensionless in the new 
variables due to the relation D x — > M 2 \3 X . Thus the action (j3J) can be written in terms 
of the new dimensionless variables as: 

2 



S — S g + Sd 



M, 



pi 



M 



We will further use a different action: 



S — S g + Sd 



d xa 



d Xy 



R 



-g 



R 



-g 



(6) 



(7) 



because the overall factor (Mpi/M) 2 is irrelevant for the classical evolution. This form 
of action is useful for further numerical analysis, since everything is dimensionless. To 
recover physical values of the Hubble parameter, energy density etc. from our numerical 
results one has to perform a straightforward inverse procedure 



M 2 Pl M 2 e, 



H^MH, t ->• M~H, 



(8) 



where dots represent the rest of the physical parameters which are relevant for our model. 
In this paper we will work with the dimensionless variables defined in Eq. (jSJ and action 
keeping the same notations for them as we had for dimensionfull variables. This will 
be assumed in every expression unless specified otherwise. 
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2.1. General formalism 

The equation of motion for the system described by the action Q can be obtained by 
the variation of S over 0: 

6 1± = Q"u 2 ^ + Q'" (V M D0) = 0, (9) 



where (') means a derivative taken with respect to □</>. It is convenient to rewrite Eq. Q 
in the following form: 

UB= -[\nQ"(B)]'(V,B)(WB), 

□0 = B. (10) 

This form will be useful for the further analysis. The system (fTOf is hyperbolic for an 
arbitrary function Q(B), thus, the Cauchy problem is well posed. In particular, this 
means that the solutions obtained for the ideal homogeneous and isotropic Friedmann 
universe are stable with respect to high frequency cosmological perturbations + . 
The energy-momentum tensor corresponding to the action (0) is 

T% = -^=^ = g, u [Q'D0 -Q + Q" (V A 0) (V A D0)] - 

- Q" [(V M 0) (V.D0) + (V,0) (V M D0)] , (11) 

It is worthwhile noting that the Energy-Momentum Tensor (EMT) does not have the 
structure of EMT of a perfect fluid. 

2.2. Dynamics in the Friedmann Universe 

As usual in inflationary cosmology, let us consider the cosmological evolution of the 
system for spatially flat Friedmann universe with the metric 

ds 2 = alt 2 - a 2 (t)c/x 2 . (12) 

The system of the equations of motion (jH)) can be written for this background in the 
following normal form 

^= - ZHB-[\nQ"{B)\B 2 , 

f-4 

ft =-*Bi + B, 
where H = a/a. The first two equations integrate out to give 

= ^-0. (14, 

In particular case of Q(B) = B 2 /2, the last expression yields that B converges to 
some constant value £>* which is not necessary equal to zero. It is precisely this 

+ It is instructive to compare this situation with that for general "fc-essence/ inflation" models where 
requirements of stability provide restrictions on the possible Lagrangians (see [II El)- 
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peculiar fact that leads to inflation in the model under consideration. In this case 
the energy-momentum tensor ([TTf has a late time asymptotic of a A term. Namely it 
is T^ v = g^Bl/2. One can expect that for a convex function Q(B) with the curvature 
Q" restricted from below the asymptotic behavior is similar. Below we will see that this 
guess is correct and will discuss the underlying reason for that. 

Up to now the system under consideration is purely kinetic because neither equation 
of motion (JHJ nor the energy-momentum tensor (fTTJ) explicitly depend on 0. Thus, in 
the case of the spatially flat Friedmann universe, the only relevant dynamical variables 
are (0, B, B). In order to complete the description of the cosmological dynamics of this 
system we need the Friedmann equation. As we have already mentioned the energy- 
momentum tensor (fTT| does not have the algebraic structure of a perfect fluid one. 
Nevertheless in the homogeneous and isotropic case of the Friedmann universe the 
only non-vanishing components are diagonal and one can formally consider the energy- 
momentum tensor as if it were of a perfect fluid. (Note, however, that the difference 
between these structures is important for the consideration of the more subtle questions 
like cosmological perturbations) . Thus we can formally define the energy density e and 
pressure p as: 

e = T° = (Q'B - Q) - Q"j>B, (15) 
p = - l -T\ = -(Q'B-Q)-Q"^B. 
Using these definitions we can write the Friedmann equation in the familiar form 



H 



2 



e 1 



V Q\B)B - Q(B) - Q"(B)<PB\ . (16) 

The equation of motion (fT3|) and the Friedmann equation (fT6|) yield the self consistent 
minimal set of equations which determines the cosmological evolution of the whole 
system - scalar field + Gravity represented by a(t). We will later refer to this 
the self- consistent case. The system is of the third order and the phase space is (0, B, B). 
As it follows from the Friedmann equation (fTH|) the energy density e is restricted to be 
nonnegative. This provides the inequality (Q'B — Q) > Q"<j)B which defines the region 
in the phase space (<f>,B,B) where the cosmological evolution is allowed (see Fig. P). 

One can see from Eq. (fT5|) that the contributions to the energy density and the 
pressure of a system can be divided into two parts: i) part (Q'B — Q) which corresponds 
to the " cosmological-term" contribution since this term enters the expressions for e and 
p with different signs; ii) part —Q"<fiB corresponding to " ultra-hard' contribution. Note, 
that the "cosmological-term" part depends only on B while "ultra-hard" part depends 
on £>, and B. The presence of the cosmological-term-like part in the expressions 
for e and p is a peculiar characteristic of the Lagrangian Q(D0). In order to obtain 
a positive " cosmological-terrri'' contribution to the energy density we should take only 
such functions into consideration for which (Q'B — Q) > 0. The last inequality is always 
true for convex functions Q(B). In this case the effective equation of state w = p/e is 
restricted to be less than one. It is worthwhile noting that in the quarters of the phase 
space where (j)B > the system is in phantom ^31 regime: w < — 1. 
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4>* < o 



Figure 1. Here we plot a spherical region around the origin of the phase space 
(</>, B, B). The colored volume represents the region where energy density e is negative. 
The surface of this volume which does not belong to the surface of the depicted sphere 
is the surface of zero energy density e = 0. Note that this forbidden regions are in the 
quarters of the {4>,B) plane where the <f>B > and consequently w < — 1. The red 
curves represent the two fixed lines given by Eq. I|17J1 , 



The system of equations (TrT^ where the Hubble parameter is expressed through 
the dynamical variables (cf>,B,B) via the Friedmann equation, possesses two kinds of 
fixed points. Namely a trivial one: (<fi, B,B) = (0,0,0) which is of no interest for our 
purposes and a nontrivial set of fixed points which form the two curves given by the 
equations 

B = 0, B = B^ L = — ^— = B * (17) 

Note that at each of these fixed points except 5* = the system is in the exact de Sitter 
inflationary regime. Further we will always assume that B* ^ 0. 

Let us now investigate the character of these fixed points. First we linearize the 
system (|T3|) around a given fixed point corresponding to a particular B*. The linearized 
system is 

dB 

— = — SH^B, 
dt 

d5B ■ , . 

nr = B ' < 18 > 
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-0.025 



Figure 2. The evolution of the dimensionless parameters <j> and B is plotted versus 
dimensionless time t. Function Q(B) is chosen to be B 2 /2. The initial values are 
= 0.1, B = 0.1, B = -0.025 



^ = - 3^*50 + [1 - 3#* 7 ] 55 + 7 £, 
where we have defined 

1 = ^mT" (19) 

and used SB = B. One can find the eigenvalues corresponding to the matrix of the 
linearized system: 

Ai = 0, A 2 , 3 = -3F*. (20) 

Because of the vanishing Ai the Lyapunov theorem is not applicable to the problem and 
we have to investigate the stability of the system in another way. Despite the failure 
of the linear analysis to prove the stability it is interesting to find a solution of the 
linearized problem. The first two equations of the system (|T8|) do not depend on the S(p 
and one can integrate them separately. The solutions are 

B = B exp(-3#*t), 

5B = 5B + ^-[l- exp (-3HJ)] , (21) 
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Figure 3. The evolution of the energy density e, pressure p and the equation of state 
w = p/e is shown on this plot. Function Q{B) is chosen to be B 2 /2. The initial values 
are cb = 0.1, B = 0.1, B = —0.025. The system quickly evolves into a nearly de Sitter 
regime. 



where B , 5B = B — are initial values for the linearized problem. From these 
solutions it is clear that the stability in the usual sense cannot be established because 
the system is degenerate and the fixed points form a curve. Thus one can expect that 
after a small deviation from a given fixed point B* the system does not return to the 
same point B* but rather to a new point B^. The natural question which arises is 
whether the new point is close to the original one. In fact from Eq. (|2TJl we obtain 
the late time asymptotic 

B„ = B o + 7^jf- (22) 

Therefore 5B* = £>« — £>* = 5B + B /3H* and for B /3H* <C 1 this deviation is 
approximately equal to the initial deviation: 5B* ~ 5B . However, as one can directly 
prove the linearized solution for 5<fi does not converge to the fixed point 5<fi = 0, if the 
point of linearization does not coincide with late time asymptotic 5**. But if we 
linearize the system around the value of the late time asymptotic B + S /3i/*, then 
the solution for 5(b is 



S(j) 



•o + tB (27-3#*) exp(-3#*t). (23) 



Thus we prove that in the linear approximation the solutions converge to the curves 
formed by the fixed points (jTT|) . The validity of the linearized solutions can be 
numerically confirmed for the appropriately chosen initial data. The numerical 
calculation reveals the following picture. The generic evolution for 0, 0, B and B 
is shown in Fig. [21 One can see from Fig. El that, indeed, the asymptotic solution is 
inflationary. The initial conditions which lead to this sort of behavior are rather general. 
The system reaches the inflationary regime relatively quickly (see Fig. and will stay 
there forever if nothing else is added to the action Q. Of course, in the real world 
there has to be a mechanism to exit from this inflationary regime. This will modify the 
evolution shown in Figs. [2] and El We will propose a natural modification and discuss 



B-Inflation 



10 



w 



-2 



-3 



-1 




1.5 



2 



2.5 



N 



Figure 4. Here we plot the equation of state w versus number of e-folds N for various 
initial data. The function Q(B) is chosen to be \B 2 . The initial date are: for the red 
curve: = 0.1, B = 0.08, B = -0.05, for the black curve: = 2, S = 0.05, B = -0.05 
for the blue one: <fi — 0.03, S = 0.5, B — 1. The equation of state w evolves to the 
nearly de Sitter value w = — 1 within approximately one e-fold. 



the details of the graceful exit in the next section, while below we will complete the 
study of the system described by (JBJ. 

2.3. The Behavior of (e,p) Trajectories 

In order to develop intuition it is useful to consider the trajectories of the system in (e,p) 
plot. A state of the system is not uniquely described by a point in the (e,p) plot. As we 
have already mentioned, to solve the full system of the equations (Friedmann equation 
(fT6j) and equation of motion (fl"3|) one has to specify the following initial conditions: 
0o, -Bo and B . In fact, we can replace these initial conditions by specifying the initial 
energy density e , the initial pressure p and the initial angle of the (e,p) - trajectory: 
tan 9 = dp/de = p/e. It is worthwhile noting that if our system had a perfect fluid like 
EMT, then the introduced angle 9 would be closely connected to the adiabatic speed 
of sound: tan# = c 2 s . From the stability requirements it would follow tan# > 0. Our 
system is free from this restriction due to the hyperbolicity of the system (JTU1) . Using the 
equation of motion (|13|) we can rewrite the formula for the initial angle in the following 
form 




uniquely define the angle 9. One can easily check that (4>,B, B) — > (e, p, 9) is one-to 
one correspondence for all finite initial values except when 0o — 0. 
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Figure 5. Here we plot the (e,p)-trajectories originating from a point where the 
universe is decelerating. The blue region cannot be reached by the system. 

Figure El represents trajectories starting with different initial angles 9 at the same 
point (eo, Po) which is located in the region (p > — e /3), where the universe is 
decelerating. Let us analyze the behavior of the trajectories. First we notice that 
in virtue of the continuity equation (|25| a trajectory can be vertical only at the point of 
the crossing of the "phantom divide" p = —e *. Namely, by the crossing from the non- 
phantom side we have 9 = Stt/2, and by the crossing from the phantom side 9 = n/2. If 
the system has a de Sitter solution as a late time asymptotic, then 9 = by approaching 
from the phantom side and 9 = n otherwise, as it follows from Eq. ifTTj) . On the other 
hand from Eq. {23$ we find that tan# = 1 can be achieved with finite values of the 
dynamical variables only if B = 0. While from the definition (fl"5|) of the energy density 
we have p = e in this case. Thus the line p = e corresponding to the ultra-hard equation 
of state is the only possible place in the (e, p) plot where tan 9 can be equal to one 
along with the values (<j), B, B) being finite. Therefore all trajectories reaching the ultra- 
hard line p = e at e ^ have the ultra-hard line as a tangent line at this point. It 
is clear that the trajectories have the angle 5/47T when touching the ultra-hard line. If 
a trajectory touches this line away from the origin (e, p) = (0, 0) then B<f) ^ and 

* The possibility of a transition through the phantom divide was discussed e.g. in Refs. [1411111 flTl] 
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the only possibility for a further evolution is to reflect in the direction of the smaller e. 
Immediately after the reflection the trajectory can only increase its angle. The angle 
is now bounded from below with 9 = 5 /An and bounded from above by 9 = 3/2tt. 
Thus, moving downwards the trajectory increases its angle and inevitably crosses the 
phantom divide. If we take into account that the angle 9 must be a continuous function, 
we find that all trajectories which cross the phantom divide and which have the initial 
angle 9 < n must reflect from the ultra-hard line. Another possibility to have tan 9 = 1 
corresponds to infinite values of B and <fi. This is a limiting case which is shown by the 
green line in Fig. For the trajectories with the initial angles large than 9 = 5/47T 
the angles are bounded from below and from above and the trajectories are evolving 
downwards. Thus, they must cross the phantom divide in any case. Considering a 
limit of the trajectories which touch the ultra-hard line we find that there must exist 
a critical trajectory which goes to the (e, p) = (0,0). For a further consideration of 
it is convenient to introduce the critical angle 9 cr = 9 cr (e , p ) which corresponds to 
the critical trajectory having the late time asymptotic: e — > 0, p — > as t — > oo. This 
critical trajectory is shown in red in Fig. [HI Then all the trajectories can be divided into 
three subclasses. 

The first subclass contains trajectories with the initial angles n/2 < 9 < 9 cr . The 
properties of these trajectories are the following: they go upwards, (if p < intersect 
the axis p = at some point), then reflect from the line p = e at the angle 57r/4. After 
that they again intersect the axis p = and go to the phantom regime crossing p = — e 
at the angle 3n/2. Finally, these trajectories approach the de Sitter regime p = — e. 

The second subclass contains trajectories with the initial angle 9 cr < 9 < 57r/4. In 
this case trajectories go directly to p = —e, approaching this line as t — > oo but never 
crossing it. This happens because these trajectories are moving downward towards the 
phantom divide with the angle bounded from above by 5/47T. A trajectory cannot have 
this limiting angle because it evolves away from the ultra-hard line and in that case the 
5/47T angle corresponds to an infinite or B. 

And, finally, the third subclass contains the trajectories with the initial angle 
57r/4 < 9 < 37r/2. These curves intersect p = — e at the angle 3n/2 going to the 
phantom regime and then approach the inflationary regime p = — e. 

Figure El shows trajectories which begin in the phantom regime (eq + po < 0) with 
different initial angles 9 . In this case the analysis is analogous to the previous case. 
The whole set of trajectories also can be divided into three subclasses as in the case of 
the "ordinary" initial conditions ( p > — e /3). The critical angle 9 cr is also defined as 
the initial angle of the trajectory which approaches e = 0,p = 0ast— > oo. The critical 
trajectory is shown in red in the Fig. El The first subclass consists of trajectories with 
initial angles 9 cr < 9 < n/2. In this case trajectories intersect p = — e at the angle 
7r/2. Then one arrives to the situation which was described by the first subclass in 
(so + Po > 0) case. The second subclass is the trajectories, which start with the angles 
7r/4 < 9 < 9 cr . The trajectories intersect the line p = —e and then one starts from the 
conditions described in the second subclass of (e + Po > 0) case. The third subclass 
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Figure 6. Here we plot the (e,p)-trajectories originating from a point where the 
universe undergoes a super-acceleration stage. The blue region cannot be reached by 
the system. 



is the set of trajectories with the initial angles — n/2 < 9 < tt/4. These are the curves 
which directly start to approach the inflationary regime p = — e. 

By the analysis performed in this subsection we confirm our result that the system 
generically approaches the de Sitter regime = — £*. The only exceptions are the 
critical trajectories (shown in red), in this case = 0. However as we have shown for 
each point (eq } Pq) there is only one solution corresponding to the critical curve. Thus 
the solutions which do not lead to the late time de Sitter regime have a zero measure. 

2.4- Stability of the system 

It is well known (see, for example, recent papers ^H]) that higher derivative theories 
possess ghost (or phantom) . The presence of a ghost degree of freedom in a dynamical 
system leads usually to various concerns about classical and quantum stability of such 
a system [IE]. While the detailed investigation of the relevancy of all these problems 
regarding our model is under way [12] we would like, perhaps, to make a few comments 
on these issues below. 
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First of all, let us briefly discuss the classical stability of the model. In the 
above analysis we have studied the dynamics governed by the homogeneous equations 
of motion only. It was shown that the solutions of these equations are stable with 
respect to the perturbations in the initial conditions. However, showing that solution 
of the homogeneous equations is stable is not enough for cosmological applications. 
Even disregarding the perturbations of the metric, one has to worry about the classical 
stability of this solution of the partial differential equations in 3 + 1 space-time. 
Regarding this issue we can note that the equations in our model are hyperbolic. 
Therefore the Cauchy problem is well posed. This means that the dependence of 
the solutions on the initial conditions is continuous and exponential instabilities are 
excluded. Thus, possible classical instabilities, if any, are not as dramatic as it is 
for systems which may flip from hyperbolic type to elliptic type in the course of the 
dynamical evolution. If one includes the perturbations of the metric, the problem 
becomes even more involved and the further analysis is needed. 

Finally, one may also worry about quantum instability due to the one graviton 
exchange between the ghost degree of freedom and the Standard Model particles. While 
during inflation the argument from [TBI i s n °t applicable, after the exit the residual 
ghost may pose such problem. There is no doubt, that weakly coupled ghost degree of 
freedom leads to the vacuum instability. However, if the ghost couples strongly to a 
regular scalar field in a particular way, which is the case in our model, the argument 
based on the evaluation of the creation rate of Standard Model particles out of vacuum 
is not reliable anymore. This issue certainly deserves further investigation. 



3. Exit Mechanism 



Every realistic inflationary scenario needs to have a natural exit mechanism. In this 
part we will modify our model in order to have an exit from the inflationary regime. Let 
us return to the action (JHJ) and add a first-order derivative term. It can be, for example, 
the first term from the Lagrangian (|TJ). Thus we shall consider the following action for 
the scalar field: 



S = d 4 xy/—g 



(26) 



M Pl J 2 \M 2 M Pl 

This choice has been widely discussed as an effective approach in string theory [17J and 
its simplicity facilitates further analysis. In a cosmological context such a term was 
considered in IE]- 

To simplify the analysis of the system we again use the dimensionless variables 
defined according to The overall factor (Mpi/M) 2 is irrelevant for the classical 
analysis and we can use the following action instead: 



S = J (fx^[K(<p)X + Qin<p)}, (27) 
where X depends on the dimensionless quantities x^ and as usual X = i(V M 0) 2 . 
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The implementation of the exit from the inflationary regime is based on the 
following idea. We assume that there are two regions of 0: one is a Q-dominated region, 
where < 1, so that the first term in (|2*Tj) does not influence the dynamics; and 

another is the region where the term K((j))X influences the dynamics significantly in 
such a way that the system exits from the inflationary regime. We will call this region 
the K + Q-region. Starting from Q-region, the system will enter the inflationary regime 
in a way very similar to that discussed in the previous section. Then during inflation 
the value of is changing as = (fio + <fi*t (where ^ according to Eq. (fTTjl ). and 
the system may come eventually to the K + Q-region. For example, if K(<p) — > for 
<p — > — oo and K(<fr) — > oo for — > +oo, then starting from inflation in the region of large 
negative <\> with > 0, our system will eventually come to the K + Q-region, where we 
expect the system to exit from the inflation. During the evolution, the contribution to 
the action which comes from the .fT-term grows during the inflationary regime whereas 
the contribution from the Q-term remains almost the same. Therefore we can roughly 
estimate the time when inflation ends: <fil\K(<j) e )\ ~ Q(B t .) where (fi e = 4>o + $*t e . The 
number of e-foldings during the inflationary regime is given by 

N~H^Z*±. (28) 

10*1 

Below we will find the necessary conditions on the function K((f>) in order to have an 
exit from the inflationary regime, and will show that a large class of functions K((f>) has 
the desired properties. 

The energy density and the pressure coming from the additional term in the action 
K((p)X give the following contribution to the energy density and the pressure of Q-part 
of Lagrangian (Eq. fT5| ): 

e 1 =p 1 = K(<p)X, (29) 

where subscript "1" reflects the fact that the term under consideration contains no higher 
order derivatives. The equations of motion ([TO become 



UB= -[\nQ"(B)]'(V,B) 2 + (Q"(B)Y 
D(p = B. (30) 
For simplicity we assume Q(B) = B 2 /2 in the further analysis, since there is no 



qualitative difference between this choice and other convex functions Q{B). The 
equation of motion (|3T)| simplifies to 

1 

2" 

□0 = B. (31) 



UB = K{<P)B + -K\<P) (V^ 



These equations of motion have the structure similar to the usual Klein-Gordon 
equation. In particular the first term on the right hand side plays the role of a mass 
term for the field B. Using this analogy we come to the conclusion that in order to avoid 
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a tachyonic instability we have to assume that K(4>) < 0. For the Friedmann universe 
the system of equations (|3TJ) takes the following form 

B + 3HB = -m 2 ((p)B - m(0)m'(0)0 2 , (32) 

(f) + 3H0 = B, (33) 

where we have introduced the effective mass m 2 (0) = —K ((/>). When the r.h.s. of 
Eq. becomes large enough then inflation ends. The value of <fi where the evolution 
changes significantly can be estimated by a comparison of the friction term in (|32|) and 
the "mass term", m 2 (<p)B. When 

H ~ m(0) (34) 

the field B becomes massive and starts to oscillate and, thus, inflation ends. This is in 
fact quite similar to the exit from inflation in "standard" potential induced scenarios. 

While the behavior of the system in the transition regime between the box- 
dominating phase and the oscillating phase is difficult to study analytically, it is possible 
to describe the system in the asymptotic oscillating regime. Our further analysis is based 
on the book [THJ chapter 4. Introducing the rescaled fields u = aB, v = a<p and the 
conformal time 77 = J dt/a, we can rewrite the system (|32|) . (|33|) in the following form: 

., (11 o!'\ v 2 dm 

u + m a ) u = "rn—rr-, 

\ a J a d(p 

v" -—v = u, (35) 
a 

where the prime denotes the derivative with respect to r/. The term on the right hand 
side in the first of these equations rapidly decays in the expanding Friedmann universe. 
Following the book [IBJ we find that if \a" /a 3 \ ~ H 2 ^> m 2 the first term in the brackets 
can be neglected as well and the approximate solution corresponds to the frozen field 
B = B*. This is exactly the inflationary solution obtained in the previous section. In 
our problem the effective mass m(<p) is growing with time, whereas the Hubble constant 
changes only slightly. As the mass becomes larger than the Hubble parameter H ~ i/* 
we can neglect the second term inside the brackets in Eq. (|35|) . The WKB solution of 
the simplified equation written in terms of B is then 



B oc a 3//2 — , sin 



m((p)dt 



v /m(0) 

This solution is valid in the case of a slowly varying mass m(0(t)), i.e. if 



(36) 



(ma)' 



m 2 a? 



$ < I. (37) 



Going back from the conformal time units to a cosmic time the above inequality becomes 

m H . . 

—z + — < 1. 38 

m z m 

Since H m already, the first term in (|38|) is also much less then one. From the 
equation it follows that 



^ J dtBa 3 oc -a~ 3/2 m" 3/2 (0) cos J m(<f>)dt 



(39) 
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Estimating the energy density we obtain 

B 2 ■ ■ m 2 6 2 1 , x 

e K T _ * T~ x (40) 

The neglected next-to-leading terms we have neglected are of the order of H(t)a~ 3 (t) 
and are oscillating. Since H{t) is a decreasing function of time and we average over fast 
oscillations these terms, indeed, quickly become subdominant. One can also estimate 
that the pressure behaves as 



p oc —a 3 cos 



m((p)dt 



(41) 



After averaging over fast oscillations one obtains that (p) =0. Thus our system has a 
dust-like behavior after the inflation is over. 

Now we are ready to formulate the conditions for K((p) in order to obtain the exit 
from inflation using the action pTjk i) K(<f>) < 0; ii) there should exist the region of (j) 
where |i^(0)| <C 1; iii) Eq. (jHll must have a solution at some e , where the value of H 
is taken as it were in the inflationary regime. The last condition means, in particular, 
that for a function K(4>) bounded below, the presence of the exit in the model depends 
on the value of B (and consequently on the value of H in f) during the inflation. Below 
we study numerically two examples of the function K(cj)) to illustrate some generalities 
of the proposed exit mechanism jj. 

The exit from the inflation which we have constructed still has to be complemented 
by a reheating mechanism, which would produce radiation after the end of inflation. 
Since the residue component has a dust-like behavior after inflation ends, the mechanism 
which converts its energy into radiation is the same as in the standard inflationary 
scenarios: we need to couple the oscillating field (in our case it can be B instead of the 
inflaton field) to some scalar fields \ through, for instance, an interaction term Bx 2 / A 
where B and \ have dimension three and one respectively, and A can be a Planck or a 
somewhat lower scale. Such an interacting term would lead to an energy density transfer 
from the massive "field" B to the radiation-like particles via, for example, the narrow 
(or broad) parametric resonance. We will discuss the details in the forthcoming paper 



3.1. Numerical examples 

First we consider the case 

K{<j>) = -e+. (42) 

We shall see below that the exit from the inflationary regime in this case is sharp 
in comparison with the next example. This may be important for generation of the 
curvature fluctuations. This will, however, be discussed in [12]. 

For the choice (JUj) the inflationary regime starts for large negative 0. There are 
two possible sets of initial conditions: <fi* < and 0* > 0. If <fi* > then 4> grows 

t) It is worthwhile noting that the discussed exit mechanism reveals some similarities with Ref. 




as (ft = cf)Q + <f>*t in the inflationary regime and the term —X becomes more and 
more important and eventually drives the system into a state where the parameter w 
oscillates as shown in Figs. [3 and |B1 The time average of these oscillations is w ~ 0, i.e. 
the system after the exit has a dust-like equation of state. The behavior of the system 
for t — > oo can be described by Eqs. l|32Jl and with m ~ e^ e,/2 , where 4> e can be 
roughly estimated using Eq. (J3IJ) . The energy and pressure are given by Eqs. ijlUJ) and 
(jlTJ) . In the case of <p* < there is no exit from the inflationary regime for the chosen 
function K((f>). 

As the second example, we take 

where a is some constant. There are several differences in the dynamics in comparison 
with the previous example. First of all, for the choice the inflationary regime is 
recovered for <fi — > +00 as well as for (j) — > —00. Another difference is that the exit 
in this example does not necessarily happen, because for sufficiently large B* the first 
term in (|27|) is not "strong" enough to drive the system out of inflation. The maximum 
value of B* that still allows an exit can be estimated using Eq. (plj) . B™ ax ~ a -1 . The 
evolution of a system for B* < B™ ax is shown in Figs. E3 and HH In addition the exit 
from inflationary stage in this case is milder than in the previous example. This is 
because the function ifO) is more widely spread than the function (jUJ. Asymptotically 
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t 




Figure 8. Energy density, pressure and equation of state for the system K{4>) — 
- exp0, Q = B 2 /2. Initial values were (f> = -20, j> = <j>* = ^/2pi, B = 10, B = 0. 



we have e — > as t — > oo, so the system arrives to the solution (|36|l with m ~ 1/a. 
4. Summary 

In this work we have introduced a new inflationary scenario where inflation is driven by 
a scalar field without a potential. We have shown that a very simple Lagrangian, which 




is a function of the second derivatives only, naturally yields an inflationary solution. 
The initial conditions which lead to inflation are generic. The inflationary solution is 
manifestly stable with respect to the cosmological perturbations. 

We have proposed a mechanism of the graceful exit from the inflationary regime 
into a dust-like state with w ~ 0. In order to generate radiation, the usual mechanism 
of reheating similar to that in the "standard" inflationary models can be applied to our 
model. The detailed analysis of reheating, cosmological perturbations and observational 
consequences will be given elsewhere [T2] . 
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